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Four Decades of Space-Borne 
Radio Sounding 



Robert F. Benson 


Abstract 

A review is given of the 38 rocket, satellite, and 
planetary pay loads dedicated to ionospheric/ magnetospheric 
radio sounding since 1961. Between 1961 and 1995, eleven 
sounding-rocket payloads from four countries evolved 
from proof-of-concept flights to sophisticated instruments. 
Some involved dual payloads, with the sounder transmitter 
on one and the sounder receiver on the other. The rocket 
sounders addressed specific space-plasma- wave questions, 
and provided improved measurements of ionospheric 
electron-density ( N e ) field-aligned irregularities (FAI). 
Four countries launched 12 ionospheric topside-sounder 
satellites between 1962 and 1994, and an ionospheric 
sounder was placed on the Mir Space Station in 1 998 . Eleven 
magnetospheric radio sounders, most of the relaxation type, 
were launched from 1977 to 2000. The relaxation sounders 
used low-power transmitters, designed to stimulate plasma 
resonances for accurate local N e determinations. The 
latest magnetospheric sounder designed for remote sensing 
incorporated long antennas and digital signal processing 
techniques to overcome the challenges posed by low 
N e values and large propagation distances. Three radio 
sounders from three countries were included on payloads 
to extraterrestrial destinations from 1990 to 2003. The 
scientific accomplishments of space-borne radio sounders 
included (1) a wealth of global N e information on the 
topside ionosphere and magnetosphere, based on vertical 
and magnetic-field-aligned N e profiles; (2) accurate in-situ 
N e values, even under low-density conditions; and (3) 
fundamental advances in our understanding ofthe excitation 
and propagation of plasma waves, which have even led to 
the prediction of a new plasma-wave mode. 

1. Introduction 

Table 1 provides a summary of the 3 8 rocket, satellite, 
and planetary payloads from eight countries (US, Canada, 
France, US SR, Norway, Japan, Sweden, and Italy) dedicated 
to ionospheric, magnetospheric, and planetary radio 


sounding, launched between 1961 and 2003. It indicates 
that the first radio sounding in space was perfonned from 
a US rocket launched on June 24, 1961, following the 
successful US rocket antenna-deployment test of June 14, 
1961. These two rockets were followed by a US rocket 
flight into a disturbed ionosphere on October 31, 1961. 

The first ionospheric sounding from a satellite was 
from the Canadian-built and US-launched Alouette 1. It 
was launched on September 29, 1962. It also established 
another first: the first satellite to be launched from the 
Western Test Range in California. It was the first satellite 
in the highly successful International Satellites for 
Ionospheric Studies (ISIS) program. The program, which 
included an additional eight countries (Australia, Finland, 
France, India, Japan, New Zealand, Norway, and the UK), 
included six satellites. Five of these satellites contained 
topside sounders. Four of the sounding satellites were 
Canadian-built and US-launched, and contained swept- 
frequency sounders (Alouette 1 and 2 and ISIS 1 and 2). 
The other two satellites in the program were built by the 
US (Explorer 20 and Explorer 31). Explorer 20 contained 
a topside sounder that operated on six fixed frequencies. 
This mode of operation proved to be so valuable for the 
investigation of N e gradients and plasma resonances that it 
was included in the design of the sounders on ISIS 1 and 2. 
Explorer 3 1 , also known as Direct Measurements Explorer- A 
(DMEA), was launched piggyback with Alouette 2. This 
dual-satellite mission was known as ISIS X. Explorer 3 1 
demonstrated the compatibility ofradio-sounder operations 
and in-situ measurements. The last satellite in the ISIS 
program, ISIS 2, was a complex observatory that included 
the first two auroral imagers in orbit. It demonstrated the 
compatibility of such imaging from a spinning satellite 
containing long sounder antennas. 

The USSR and Japan launched an additional seven 
topside sounding satellites (Cosmos 381, ISS A and B, 
Intercosmos 19, EXOS C, Cosmos 1809, CORONAS 1) 
between 1970 and 1994. Cosmos 381 followed the launch 
of Cosmos 3 1 8 in 1 969 that failed due to a rocket problem. 
The Cosmos 381 sounder, which was a pulsed sounder 
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Table 1. A summary of space-borne radio sounders*. 


Date 

Country 

s/c 

T 

Freq 

Dp (m) 

P(W) 

Sounder Experiment and Comments 

Sample Results and Reference 

06/14/61 

US 

r 



75 

no tx 

Antenna deployment test from spinning s/c 

Successful demonstration [ 1 ] 

06/24/61 

US 

r 

s 

f 

32 

3 

4.07 & 5.97 MHz, 22 pps, 1,060 km apogee 

1 st topside sounding, quiet [2] 

10/31/61 

US 

r 

s 

f 

32 

3 

4.07 MHz, 22 pps, 1,070 km apogee 

t/s sounding of spread F [3] 

09/29/62 

Canada 

s 

s 

s 

23,46 

100 

Alouette 1, (0.5-12) 63 pps, 1,000 km polar 

1 st Global N p (h ) t/s profiles [4] 

08/25/64 

US 

s 

s 

f 

19, 37 

45 

Exp. XX, 6 freq (1.5 to 7.22), 3 dp at 60° 

N e & res at 800 km polar [5, 6] 

11/29/65 

Canada 

s 

s 

s 

23, 73 

300 

Al. 2, (0.1-14.5), 30 pps, 500-3,000 polar 

Compare in-situ techniques [4] 

01/30/69 

Canada 

s 

s 

f/s 

19, 73 

400 

ISIS 1 , (0. 1 -20), 30 pps, 550-3,500 polar 

More high-alt sound., AKR [4] 

10/06/70 

France 

r 

r 

s 

4 

low 

EIDI 1, 1-11 MHz in 0.4 s, 500 pps 

3 f„ res: 2 waves 1 kHz apart [7] 

10/22/70 

France 

r 

r 

s 

4 

low 

EIDI 2, 2 nd of 3 rocket relax sounders 

Results given for EIDI 1 & 3 [8] 

12/20/70 

USSR 

s 

s 

f 

18 

100 

Cosmos 381, 985-1023 km, 74° inclination 

N e irreg. related to obs. UV [9, 10] 

04/01/71 

Canada 

s 

s 

f/s 

19, 73 

400 

ISIS 2, (0.1-20), 45 pps, 1,400 km polar 

Compatible/w imagers, etc. [4] 

06/25/71 

US 

r 

r 

s 

10 

60 

0.5-4. 5 MHz in ~ 4 s, 50 pps, waveform tel 

f T res freq variation obs. [11] 

01/13/72 

Norway 

r 

r 

s 

5 

low 

Mother-Daughter, 1.1-9 MHz, 192 km apog 

Res identified as cone res [12] 

02/06/73 

Norway 

r 

r 

s 

8 

low 

Mother-Daughter, 0.8-5 MHz, 260 km apog 

Cone & other res observed [13] 

02/17/73 

France 

r 

r 

s 

4,8 

2 

EIDI 3, 0.44-5.52 MHz in 4 s, 10 ms listen 

N p irreg from f N res [14] 

02/29/76 

Japan 

s 

s 

s 

11,37 

n/a 

Ionosphere Sounding Satellite (ISS)-A 

One month lifetime 

04/20/77 

France 

s 

r 

s 

42 

low 

GEOS 1, 0.3-77 kHz in 256 300 Hz steps 

Magnetospheric N p , 5-7 R c , [15] 

10/22/77 

France 

s 

r 

s 

215 

low 

1SEE 1 (mother) w/sounder & 2 (daughter) 

Magnetospheric N e [16] 

02/16/78 

Japan 

s 

s 

s 

11,37 

150 

ISS B, (0.5-14.8), 9 & 64 pps, 972-1220 km 

Global maps of foF, [17] 

07/14/78 

France 

s 

r 

s 

42 

low 

GEOS 2, 0.3-77 kHz in 256 300 Hz steps 

Plasmapause N e features [18] 

09/16/78 

Japan 

s 

r 

s 

73,102 

300 

EXOS B (Jikiken), (.02-3), 227-30051 km 

Plasmapause N e irreg. [19] 

02/27/79 

USSR 

s 

s 

s 

15, 50 

300 

Intercosmos 19, 0.3-15.9 MHz in 5.8 s 

Sounder-accel. electrons [20] 

02/14/84 

Japan 

s 

s 

s 

40 

300 

EXOS C (Ohzora), (0.1-16.0), 354-865 km 

Wave-mode conversion [21] 

02/22/86 

Fran/Swed 

s 

r 

s 

80 

low 

V4H on Viking, 8-500 kHz in 3 freq. bands 

N e , plasmasphere & polar cap [22] 

12/18/86 

USSR 

s 

s 

s 

15,50 

300 

Cosmos 1809, (0.3-15.95), 58.6 pps 

Freq. of sound-ace electrons [23] 

01/30/89 

Canada 

r 

s 

s 

4.4, 5.5 

2 

OEDIPUS A, 0-5 MHz bi-static sounder 

Guided Z-mode propagation [24] 

02/22/89 

Japan 

s 

s 

s 

60 

600 

EXOS D (Akebono) 0.3-11.0, 10500 km ap. 

Magnetosph. ducted echoes [25] 

10/06/90 

US/France 

p 

r 

s 

72 

low 

URAP on Ulysses, (0-50 kHz) 

N p in lo plasma torus [26-28] 

03/03/94 

USSR 

s 

s 

s 

15 

n/a 

CORONAS 1, 500 km, 82.5° inclination 

Upper-hybrid emission band [29] 

11/07/95 

Canada 

r 

s 

s 

4.4, 5.5 

10 

OEDIPUS C, 0.1-8 MHz bi-static sounder 

In-situ Faraday-rotation [30] 

10/15/97 

US/France 

p 

r 

s 

20 

low 

RPWS on Cassini, 3.6-115.2 kHz in 90 steps 

N p in vicinity of Saturn [3 1 ] 

O 

oo 

VO 

00 

USSR 

s 

s 

s 

n/a 

n/a 

Orbital Complex (OC) on Mir Space Station 

Sounding below F peak [32] 

03/25/00 

US 

s 

s 

s 

500, 20 

10 

IMAGE/RPI, (.003-3 MHz), tuned tx ant. 

2D F-AN structures [33, 34] 

07&08/00 

France 

s 

r 

s 

88 

low 

WHISPER on 4 Cluster s/c, 2-80 kHz 

Multipoint magnetospheric N e [35] 

06/02/03 

Italy/US 

p 

s 

s 

40 

5 

MARSIS on Mars Express 

Radio sounding of Mars iono [36] 


*The date is the month/d ay /year of launch; s/c indicates that the spacecraft was either a rocket (r), satellite (s), or space probe 
(including a satellite around another planet) (p); T indicates the type of sounding, i.e., either a low-power relaxation sounder (r) 
for resonances or long-range sounding (s); Freq indicates that the sounding was at one or more fixed frequencies (f or was a 
swept-frequency sounder (s); Dp (m) indicates the tip-to-tip length of the sounder dipole antenna(s) in meters; P (W) indicates the 
estimated peak radiated pulse power in watts; n/a indicates that the information could not be located by the author. 


using 20 fixed frequencies, was noteworthy for at least 
two reasons. First, the sounder design team was awarded 
the silver medal of the Moscow Exhibition of National 
Achievements. Second, one member of the design team 
(A. I. Galkin) was the father of I. A. Galkin, who is a key 
member of the sounder on the IMAGE satellite known as 
the Radio Plasma Imager (RPI), to be discussed below. In 
1998, the USSR placed a radio sounder on the Mir Space 
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Station. It was the first ionospheric sounder to be placed in 
orbit near the altitude of the N e maximum, h max . 

From 1970 to 1995, eight ionospheric-sounding 
rockets built by four countries (France, US, Norway, and 
Canada), including four dual payloads (two each from 
Norway and Canada) were launched to investigate wave 
propagation and plasma resonances. The two French 
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rocket sounders launched in October 1970 were the first to 
telemeter the sounder-receiver waveform to the ground for 
plasma-resonance investigations. The US rocket sounder 
launched in June 1971 also transmitted the sounder-receiver 
waveform to the ground, it was the first “smart” sounder, 
in that the mode of operation was changed when strong 
plasma resonances were encountered. 

From 1977 to 2000, eleven magnetospheric sounders 
(on the GEOS 1, ISEE 1, GEOS 2, EXOS B, Viking, 
EXOS D, IMAGE, and four Cluster satellites) from four 
countries (France, Japan, Sweden, and the US) were placed 
into orbit. The first was of the relaxation type, i.e., low- 
power transmissions designed to stimulate local plasma 
resonances. It operated from the European GEOS 1. Since 
the intended geostationary orbit was not achieved, most 
local electron-density ( N e ) measurements, as deduced 
from the sounder-stimulated plasma resonances, were 
made from radial distances between five and seven Earth 
radii ( R [: ). The last four were also relaxation sounders, 
designated as WHISPER (Waves of High frequency and 
Sounder for Probing Electron density by Relaxation); an 
identical sounder was placed on each Cluster satellite. The 
first magnetospheric sounding (out to a radial distance 
of nearly 2 R E ) to stimulate the full spectrum of plasma 
resonances, including a class of resonances known as the 
Dn resonances that were prominently stimulated in the 
ionosphere by Alouette 2 and ISIS 1 at radial distances 
out to about 1.5 R E , was from the Japanese EXOS D 
(Akebono). The first long-range magnetospheric sounding 
(out to a radial distance of 8 R E ), using a programmable 
instrument with digital signal processing, was performed 
by the US sounder RPI on the IMAGE satellite. It was 
capable of low-frequency operation because of 500-m 
tip-to-tip dipole antennas in the spin plane. As in the case 
of the ISIS-2 observatory, IMAGE/RPI demonstrated the 
compatibility of remote imaging from a spinning satellite 
containing long sounder antennas. 

Three extraterrestrial sounders, based on cooperative 
agreements among the US, F ranee, and Italy, were launched 
between 1990 and 2003. Two were of the relaxation type: 
the Unified Radio And Plasma- wave instrument (URAP) on 
the Ulysses spacecraft to Jupiter (encountering the Io plasma 
torus in February 1992), and the Radio and Plasma Wave 
Science (RPWS) investigation on the Cassini spacecraft 
to Saturn. The Mars Advanced Radar for Subsurface and 
Ionospheric Sounding (MARSIS) instrument performed 
the first planetary ionospheric sounding at Mars in June 
2005 from the ltalian/US Mars Express. 

While the above discussion highlights the many firsts 
associated with space -borne sounders, it is important to note 
that there was a fundamental first associated with ground- 
based sounders, commonly referred to as ionosondes. As 
pointed out by Warren Flock [37], 

The ionosonde, developed in 1925-1926 and used for 

monitoring the ionosphere since, was the first practical 
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radar system of significance to be put into service. 
Operating in the HF frequency range and developed 
some years prior to surveillance radars, before the term 
radar (for radio detection and ranging) was introduced, 
the ionosonde may not be thought of by some as a radar, 
but it clearly is a special-purpose radar system. 

Descriptions of the operation of swept-frequency 
ionosondes, and the interpretations of their data records, 
called ionograms, have been given in numerous books and 
publications [38-45]. Modern ionosondes, in use worldwide, 
incorporate both enhanced remote-sensing capabilities and 
automatic-analysis techniques [46-51]. 

The material in Table 1 was based on a comprehensive 
review of the ISIS program [52], a more-recent review that 
includes sounders other than those from the ISIS program 
[53], and information collected by the author with the aid 
of colleagues. It includes representative references from 
the various space -borne sounders: see [54] for a thorough 
review with references covering the first two decades of 
Alouette-1, Explorer-20, Alouette-2, and Explorer-31 
operations. The related discussion in the text is based on the 
author’s involvement with radio-sounding research using 
data from sounders carried on the Alouette 1 and 2, ISIS 1 
and 2, OEDIPUS C, Ulysses, and IMAGE spacecraft. The 
objectives, history, and principal achievements of the ISIS 
program were discussed in a special issue ofthe Proceedings 
of the IEEE on topside sounding and the ionosphere [55]. 
There has been a resurgence of interest in the analysis 
of the Alouette-2, ISIS-1 and ISIS-2 topside sounder 
data, following a data-rescue effort that produced digital 
ionograms from a subset of the original analog telemetry 
tapes - many that were not previously processed into 
35-mm film ionograms [56]. OEDIPUS C (Observations 
of Electric-field Distributions in the Ionospheric Plasma -a 
Unique Strategy) established a new record length for the 
longest tether in space at the time (approximately 1.2 km), 
breaking the record set by its predecessor, OEDIPUS A. 
The tether was cut shortly after the dual rocket payload 
achieved apogee, in order to provide the setting for unique 
bistatic radio-sounding experiments. The addition of an 
active sounding mode to the URAP instrument on the 
interplanetary Ulysses probe was the result of leadership 
by the late R. G. Stone, as PI, after instrument selection. 
IMAGE/RPI employed the longest antenna elements ever 
deployed on a spinning spacecraft, namely, orthogonal 
500 m tip-to-tip dipoles in the spin plane, and a 20 m tip- 
to-tip dipole along the spin axis. The former were later 
shortened (presumably by micrometeorite impacts), but 
did not prevent successful active soundings (see, e.g., 
Section 1.2 of [57]). 

This paper is not intended to be a thorough review 
of space-borne radio sounding, as approximately 1,000 
scientific papers have been written based on these sounders. 
Rather, the plan is to discuss space-home radio sounding as 
the gold standard for in-situ and remote electron-density ( 
N e ) determinations in Section 2. Fundamental plasma 
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Figure 1. (a) An Alouette-2 mid- 
latitude ionogram, in negative 
format with signal reception in 
white on a black background, 
showing Z, O and X traces, 
plasma resonances, and a natural 
noise band (presumably in the 
Z mode): see the discussion in 
the text, (b) The group veloci- 
ties corresponding to the plasma 
conditions at the satellite, (c) The 
calculated N e values from each 
of the traces. (Figure adapted 
from [59, 60].) 
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processes, and gradients in N e and the magnetic field, B, are 
discussed in Section 3. Enabling instrumental innovations 
are covered in Section 4. The importance of the antenna’s 
orientation (relative to B) and of special plasma conditions 
for sounder-stimulated plasma phenomena are discussed in 
Sections 5 and 6, respectively. Nearly lossless propagation 
within N e wave ducts are covered in Section 7. Section 8 
provides a brief summary. 

2. Space-Borne Radio Sounding 
as the Gold Standard for N e 
Determinations 

The original motivation for placing radio sounders in 
space was to obtain routine ionospheric N e information 
over large geographic regions at altitudes above h max , i.e., 
in the topside ionosphere: hence, the name ionospheric 
topside sounders. They are complementary to ground-based 
ionosondes, in that they provide orbit-plane topside N e 
contours over a short time interval above a ground-based 
station where an ionosonde provided 24-hour coverage 
of the N e profile below h max . (Note: In each case, there 
are problems associated with the determination of h max . 
Ionosonde measurements from the ground typically yield 
values that are about 10 km too low, and topside-sounder 
measurements from a satellite typically yield values that are 
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about 10 km too high [58].) Since many of these topside- 
sounder satellites were in high-inclination orbits, the 
topside orbit-plane N e contours corresponded to latitudinal 
contours. With an orbital period of about 90 minutes, these 
contours can be considered to be snapshots of the topside 
ionosphere covering about 40° of latitude in 10 minutes. 

In addition to the differences in spatial and temporal 
coverage between ionosondes and topside sounders, there 
are fundamental differences in the received signals. These 
signals are recorded on data records called ionograms, 
which display the apparent range to the reflection point as 
a function of sounding frequency. The apparent range, also 
called the virtual range, is defined as ct/2 , where c is the 
vacuum speed of light, and t is the roundtrip travel time. 
The apparent range is greater than the true range because, 
in the frequency range of main interest, the waves travel 
with a group velocity, v g , considerably less than c. An 
example of a topside ionogram is presented in Figure la. 
The values of v g for each of the principal ionospheric 
reflection traces, corresponding to the plasma conditions at 
the satellite as determined from the resonances and cutoffs 
(see next paragraph) on the ionogram of Figure la, are 
presented in Figure lb. 

Topside ionograms, such as the one in Figure la, 
differ from ground-based ionograms in that: 
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1. The apparent ranges of the echoes are much greater 
(they are typically within 1,000 km on ground-based 
ionograms [42]). 

2. The echo traces extend to zero range (they typically 
do not appear below 100 km range on ground-based 
ionograms [42]). 

3 . The Z-mode wave is directly observed (it is only detected 
after a coupling to the O-mode wave on ground-based 
ionograms [39]). 

4. Prominent plasma resonances, appearing as stalactites 
hanging from the zero apparent-range baseline, are 
present (these signals are mainly attributed to the 
reception of short-range electrostatic-mode echoes, and 
thus are not observed by ground-based sounders, where 
the antenna is in the neutral atmosphere [61]). 

5. The ionosphere often generates noise bands that are 
received over the entire listening range, such as the 
one with a low-frequency cutoff coinciding with fzS 
in Figure la. 

6. The ionosphere provides shielding of manmade 
interference in the frequency region below the 
ionospheric penetration frequency, fxF2 [62], allowing 
low powers to be used in topside sounding (a few hundred 
watts, rather than several kW as used in ground-based 
sounding in the 1960s), even though the propagation 
distances are much greater. 

7. In the frequency region above fxF2, surface reflections 
are often received that enable information to be deduced 
about the entire (topside and bottomside) N e (h ) profile 
[63], as well as terrain information [64]. At Mars, the 
sounder has been used as a ground-penetrating radar in 
this frequency region [65]. 

The labels at the top of the ionograms in F igure 1 a use 
notation common in ionospheric publications. The labels 
correspond to (with commonly used magnetospheric physics 
notations in parentheses) the electron cyclotron frequency, 
f H ( f ce ) and harmonics nf H ( nf ce ); the electron plasma 
frequency, f N (f pe )', the upper-hybrid frequency, f T ( 
f uh ); the cutoff frequency of the free-space extraordinary 
(X) mode, f x S (the “S” is used to designate the cutoff 
frequency at the satellite); and the cutoff frequency, f z S 
, of the slow branch of the X mode, called the Z mode. 
The cutoff frequency of the free-space ordinary (O) mode 
corresponds to f N . These frequencies have been defined 
in many books and review articles, e.g., see Sections 2 of 
the reviews [66] and [67]. They are defined below in terms 
of how they can be used to determine |fi| and N e : 

|B| [nT] s 35.7 (7?/^[kHz])/«, where n = 1 , 2 , 3 ,... ( 1 ) 


N e [ C m~ 3 ]^(f N [ kHz]) 2 /80.6, (2) 
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f N =(fT~l h) V2 , ( 3 ) 

/n = \_fx (fx - /h )] > ( 4 ) 

fx = [/z (/z + ///)] ■ ( 5 ) 

Equations (4) and (5) lead to the useful expression 

fH = fx ~ fz ■ ( 6 ) 


The quantity f H (and hence B from Equation (1)) 
can often be determined to an accuracy of a few tenths of 
a percent, since there are typically more nf H resonances 
present on an ionogram than in the example of Figure la, 
and the individual frequencies can typically be determined to 
an accuracy of a few percent or better (see, e.g. , Appendix A 
of [56]). Once f H is accurately determined, Equations (3) 
to (5) provide three independent measurements of f N (and 
hence, N e from Equation (2)), based on the observed values 
of f T , fx , and f z . Thus, together with the observed 
value for the f N resonance and Equation (2), there are 
four independent determinations of N e available from an 
ionogram such as the one shown in Figure la. Even when 
X and Z traces are not present, other resonances are often 
present to enable a spectral classification so as to determine 

_'3 

N e values to a precision of ~1%, even when N e ~ 1 cm 
(see Section 3, and Section 3 of [68]). 

The above redundancy in the determination of f N , the 
typical accuracy of a few percent or better in determining 
the frequencies, and the minimal effect of spacecraft/plasma 
interactions on the determined values (see Section 3), leads 
to the concept of the space -borne sounder as the gold 
standard for the accurate determination of N e . Indeed, 
the comparison of different DMEA probe techniques used 
the nearby Alouette-2 topside-sounder as the N e reference 
level [69, 70]. Also, a comparison of data acquired from 
the DMEA and Alouette-2 cylindrical electrostatic probes 
indicated “no inherent incompatibility” between sounder 
and probe operations [71]. More recently, interleaved active 
and passive Radio Plasma Imager observations were used to 
confirm that accurate magnetospheric N e measurements - 
i.e., within a fewpercent-canbe deduced from the spectral 
features of the upper-hybrid emission band, which is often 
present on magnetospheric dynamic spectrograms [72]. 
On the other hand, in the ionosphere earlier investigations 
based on comparing passive wave observations and active 
sounder operations have shown that intense wave emissions 
seldom correspond to the upper-hybrid band, particularly 
to emissions at f T (or f uh ) [29, 73] 

However, the main objective of an ionospheric topside 
sounder is not to obtain accurate in-situ N e values, but 
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to obtain complete vertical N e profiles, N e [h ), from 
the satellite altitude down to h max . This information is 
obtained by inverting the following integral expression 
for the apparent range as a function of frequency, P'(f ) 

, in order to obtain the electron-density distribution as a 
function of true altitude, i.e., N e (h ) : 

P'(f) = \n'dh, (7) 

where «' = «'{/V e (/;),/, |B|(/z),^(/z)} is the group 
refractive index, <f> is the angle between B and the direction 
of propagation, and the integration is from the satellite 
altitude down to the height of reflection for frequency f 
Equation (7) and its solution have been discussed in many 
texts and papers. The approach commonly used in the ISIS 
program is that of Jackson [60, 74, 75]. The approach is 
to partition the total integration path in Equation (7) into 
a number of laminations equal to the number of points 
selected on the X-mode reflection trace. The values of |fi| 
and <f> at the top of each lamination interval, including 
the first one starting at the satellite altitude, are obtained 
from a magnetic-field model, the variation of |fi| within a 
lamination is based on the inverse cube law. 

While the observed nf H plasma resonances are 
not used to improve on the model’s f H value, other 
resonances are used with this model value to obtain a better 
determination of f x S in situations when it is difficult to 
accurately measure (see Section 3 of Appendix B of [74]). 
In some cases, particularly when N e at the satellite altitude 
is very low, this iteration technique fails to converge. A 
modified iteration method has been shown to converge 
in such cases [76]. Refinements to true-height inversion 
techniques have continued (see, e.g., [77, 78]), andprograms 
have been developed for the automatic processing of 
topside digital ionograms [79-82], which have been recently 
produced for Alouette 2, ISIS 1, and ISIS 2 [56]. 


Since the expressions for n in Equation (7) are 
different for the O and X modes, and since the X mode 
has two branches, corresponding to a fast wave (termed the 
X mode) and a slow wave (termed the Z mode), the three 
ionospheric reflection traces, i.e., the Z, O, and X traces in 
Figure la, yield three independent ways to determine the 
N e ( h ) profile. The resultsofapplyingtheinversionanalysis 
to these three reflection traces in Figure la are shown in 
Figure lc.Thegoodagreementamongthecalculated N e (/?) 
profile values based on the observed P'(f ) values from 
the Z, O, and X reflection traces provides confidence that 
these traces were properly identified, and that the two main 
assumptions used in the analysis were satisfied, namely, that 
the propagation was vertical and that the ionosphere was 
horizontally stratified. Note that the points corresponding 
to the O and X modes extended all the way down to near 
h max , i.e., the altitude of maximum N e near 200 km, 
whereas the points corresponding to the Z mode extended 
only down to about 500 km. This limitation was caused by 
the restriction of the Z mode to frequencies less than f z I 
(called “Z infinity”) in Figure lb.lnFigure la, f z I appears 
as a narrowband signal just beyond 0.9 MHz hanging from 
the top ofthe ionogram, and it coincides with the asymptotic 
Z trace at the bottom of the ionogram. f z I varies between 
f N and f T , depending on the angle of propagation. In spite 
of this propagation limit for the Z-mode echoes, IMAGE/ 
RPI Z-mode echoes have been used to determine accurate 
N e profiles along B for thousands of kilometers above the 
IMAGE satellite, i.e., into regions of decreasing N e , when 
IMAGE was in the inner plasmasphere and at moderate to 
low altitudes over the poles [83]. These Z-mode echoes 
are possible because the Z-mode cutoff frequency, f z , in 
Equation (5) can have a local minimum leading to a trapping 
region for Z-mode waves [84]. An inversion routine was 
developed to convert the IMAGE/RPI Z-mode echoes into 
field-aligned N e profiles in the vicinity of the spacecraft 
[ 83 ] . It was shown that these profiles were in good agreement 
with model values based on a different inversion procedure 
[85] used on IMAGE/RPI X-mode echoes . 



Figure 2. ISIS-1 sounder-derived N e 
contours (TO 2 . 3xl0 2 , 10 ^ , etc.) 
and field-aligned projections through 
symbols indicating satellite position and 
auroral kilometric radiation intensity, 
with upper bold portions indicating the 
depth beneath the satellite of the auroral 
kilometric radiation source region (with 
estimated error limits), as determined 
from the spectral range of the auroral 
kilometric radiation (0959-1008 UT on 
19 August 1969; adapted from [92]). 
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In addition to Z-, 0-, and X-mode echoes, space- 
borne sounders have also detected echoes from waves 
propagating in the whistler ( W) mode. Whistler-mode echoes 
were detected on high-latitude, early-morning Alouette-2 
ionograms, and were attributed to ground reflections and 
a mode-coupling process in the lower ionosphere [86]. 
IMAGE/RPI commonly observed W- and Z-mode echoes 
when IMAGE was in the innerplasmasphere and at moderate 
to low altitudes over the polar regions [87]. Two types 
of W-mode echoes were observed by the Radio Plasma 
Imager: discrete and diffuse. The former were attributed 
to reflections from the bottom side of the ionosphere; the 
latter were attributed to scattering from field-aligned N e 
irregularities (FAI) located within 2,000 km Earthward 
of the satellite, and near the magnetic field line passing 
through the satellite. These echoes were consistent with 
cross-B field-aligned irregularity scale sizes ranging from 
10 m to 100 km. They support previous evidence [88] for 
the common occurrence of field-aligned irregularities in 
the topside ionosphere. Both types of W-mode echoes, and 
diffuse Z-mode echoes, have great potential for local and 
remote N e diagnostics [87]. 

The results displayed in Figure lc confirm the 
validity of the inversion procedure on a given ionogram 
with good traces. Furthermore, high-latitude conjunction 
studies between different topside-sounding satellites show 
good agreement among the N e (h ) (field-aligned) profiles 
obtained from the different sounders [89]. However, 
accumulated evidence indicates that the topside Alouette-1 
N e ( h ) profiles are too low, i.e., that they should be shifted 
upward parallel to the altitude axis [58]. This evidence 
included comparing simultaneous Alouette-1 topside 
profiles, ground-based incoherent scatter radar, and rocket 
measurements. It indicated that the profiles were about 
20 km too low at altitudes below 600 km. Detailed topside - 
sounder/bottom-side-ionosonde comparisons, involving 
both Alouette-1 and Alouette-2 profiles, indicated that 
the topside altitude error was only a few percent of the 
propagation path. The differences were too large to be 
attributed to the vertical propagation assumption (rather 
than performing actual ray tracing). When the propagation 
path was short, i.e., corresponding to low satellite altitudes, 
a systematic error was detected in the range markers. An 
inspection of the ground reflection traces observed on 
Alouette-2 ionograms indicated that the best N e (h ) profile 
results were obtained by assuming this systematic error 
to be 30 km, i.e., by reducing the virtual range by 30 km 
[58]. This range-marker problem led to the addition of a 
calibration range marker at 1667 km on the ISIS-2 topside 
sounder [90]. Later work, involving magnetic-field-aligned 
conj unction comparisons between low-altitude DE-2 in-situ 
Langmuir-probe measurements and ISIS-1 and -2 topside 
N e ( h ) profiles, indicated that the sounder-derived N e 
values agreed with the probe measurements to within 30% 
near h max on three of four comparisons. The agreement 
for the fourth comparison, in a region of N e irregularities, 
was within 60% [91]. 
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As stated by Jackson [58] concerning the errors 
discussed above, they “are usually too small to detract 
significantly from the general usefulness of topside 
ionograms.” Many investigations have indicated this 
usefulness of the topside N e (h ) profiles. For example, 
Alouette-2, ISIS-1, and ISIS-2 topside N e ( h ) profiles were 
used to reveal a tongue of F-region ionization (attributed to 
anti-sunward drift caused by magnetospheric convection) 
extending from the dayside across the polar cap during 
a magnetically-quiet day [89]. Topside ISIS-1 N e (h ) 
profiles were combined to produce orbit-plane contours 
that together with simultaneous passive observations of 
the auroral kilometric radiation ( AKR) frequency spectrum 
by the sounder receiver, determined the altitude extent of 
the auroral kilometric radiation source region below the 
satellite [92] (see Figure 2). Such orbit-plane N e contours 
were compared with plasmapause magnetic-field lines 
(determined near the same time period by other satellites) 
to provide examples where the plasmapause projection 
occurred near the sharp low-latitude boundary of the 
nighttime ionospheric N e main (mid-latitude) trough 
at high altitudes, but near the N e minimum at lower 
altitudes. These contours were derived from ISIS-2 N e (h ) 
hand-scaled profiles obtained from the National Space 
Science Data Center (NSSDC) ftp site. They may actually 
correspond to magnetic-field-aligned profiles, making the 
stated altitudes on the contours slightly low (see Figure 4 
of [93]), but not impacting the conclusions of the study. 

Large satellite topside sounding databases have 
enabled a number of diverse investigations. Alouette-1 and 
2 hand-scaled N e (h ) profiles have been used to create 
global-averaged N e ( h ) distributions [94]. Computer- 
scaled ISIS-2 digital ionograms produced N e ( h ) profiles 
that have been used to improve the International Reference 
Ionosphere (IRI) [95], and that have been used with IMAGE/ 
RP1 high-latitude magnetic-field-aligned magnetospheric 
N e profiles to model the topside ionospheric F region and 
plasmasphere [96]. The International Reference Ionosphere 
is in particular need of improvement in the high-latitude 
topside ionosphere, as illustrated in Figure 10b of [56]. 
Global averaged F-region peak critical frequency (foF2) 
distributions obtained from the IS S-B topside sounder have 
been used to evaluate global foF2 models [97]. Regional 
averaged foF2 and h max distributions obtained from the 
Intercosmos- 19 topside sounder have been used to detect 
ionospheric signatures of tectonic activity [98], and for 
event studies, such as identifying ionospheric precursor 
signatures of large earthquakes, as reviewed in [53]. These 
are only a few examples of the scientific results that have 
come from the remote-sensing capabilities of space-borne 
radio sounders, particularly in the ability to derive topside 
N e ( h ) profiles. More examples can be found in a detailed 
review [54] of the results from the first two decades of the 
early satellites in the ISIS program. 

The next sections will concentrate on some of the 
serendipitous scientific results that have resulted from 
these sounders, i.e., results not directly related to the main 
objective of obtaining vertical N e profiles. 
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Local es mode echoes 



Figure 3. A schematic illustration of ray paths leading 
to echoes (at times t j and t-> ) of sounder-generated 
(at time t$) electrostatic (solid curves) and electro- 
magnetic (dashed lines) waves [103], 

3. Fundamental Plasma 
Processes, and Gradients 
in N e and B 

Spectacular signal returns on topside ionograms can 
oftenbe explained using ray tracing, and special properties of 
the refractive index under certain plasma conditions, rather 
than invoking unusual N e distributions. A classic example 
is the observation of a kink in the Z reflection trace that has 
been explained in terms of the peculiar properties of the 
Z-mode group refractive index when the wave frequency 
varies near f [99]. Another is the concept of oblique 
echoes of sounder-generated electrostatic waves, due to the 
sensitivity of the refractive index to gradients in N e and 
B near / , f uh , and nf ce , to explain plasma resonances 
observed at these frequencies [100-102] (see Figure 3). 
Since these waves travel hundreds to thousands of meters 
into the medium, geophysical parameters deduced from 
the resonant frequencies using Equations (1) to (5) are not 
significantly affected by spacecraft/plasma interactions. 

Note that the ray paths for the nf ce resonances in 
Figure 3 exclude the case for n = 1 . The parabolic nf ce 
ray paths illustrated in Figure 3 correspond to electrostatic 
waves, with wave vectors, k, nearly perpendicular to B. The 
group velocities, \ s , for these waves are approximately 
parallel to B when the wave frequency, f, is near nf ce . 
Observational and theoretical investigations support the 
concept illustrated in Figure 3 for the resonances with 
n > 1 . Additional effects may play a role in the case of 
n = 2 , and the ray paths of the higher-order resonances in 
the topside ionosphere may be very short (see [104] and 
references therein). The nf ce resonances stimulated by 
magnetospheric sounders are observed for much longer 
time durations [15]. The oblique-echo mechanism has yet to 
be applied to explain these observations. The nf ce plasma 
wave dispersion curves with k approximately perpendicular 
to B are quite different for the case n = 1 , however, and an 
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interpretation of the observed features of the f ce resonance 
has been sought for decades. A recent explanation for this 
resonance involves a solution to the hot-plasma dispersion 
equation with /near f ce when k is approximately parallel to 
B [ 1 05] . It can explain the fringe pattern observed on the f ce 
resonance during fixed-frequency operations in terms of the 
beating of two waves radiated from the plane containing the 
antenna and B. The waves in this interpretation are highly 
damped, which may explain the observed strong dependence 
of the duration of the f ce resonance with altitude (longer 
duration under low-density conditions at high altitudes) in 
the Alouette-2 topside-sounder data [106]. 

There are two prominent sequences of resonances 
oftenobservedbetweenthe nf ce resonances that are not well 
explained by the oblique-echo mechanism. One is observed 
in the frequency region above f uh [107], and the other is 
in the frequency region below f pe [108], following the 
earlier discovery of the first member of the sequence [ 1 09]. 
The former are known as the Qn resonances, and have been 
explained in terms of a near matching of slowly -propagating 
electrostatic waves (known as Bernstein-mode waves) 
and the satellite velocity, [110]. The latter are known 
as the Dn resonances, and have motivated five different 
theoretical interpretations (see the review in Section 1 
of [ 1 1 1 ]). One proposed that the Dn resonances provided 
evidence for a new wave mode involving eigenmodes of 
cylindrical-electromagnetic plasma oscillations [112, 113] 
based on earlier theoretical predictions, involving force- 
free electromagnetic field (FFEMF) configurations, and 
that such oscillations would lead to a quantized emission 
spectrum with the frequencies of the different elements 
proportional to sfn [114, 115]. The cylindrical structures 
in this model correspond to the commonly observed field- 
aligned irregularity observed in space, not to the sheath 
surrounding the sounder antenna. The wave nature of 
these cylindrical oscillations was predicted as force-free 
electromagnetic field waves and solitons [116] that were later 
observed in a laboratory plasma [117]. The Dn resonances 
have prominent subsidiary resonances [108] that have been 
interpreted differently [118-120]. In the magnetosphere, 
the frequencies of the Qn resonances have been observed 
to deviate from the values expected for a Maxwellian 
electron velocity distribution [15, 68, 121, 122]. These 
observations have been investigated theoretically in two 
ways: (1) in terms of two Maxwellian components [121]; 
and (2) as a kappa distribution, where observations were 
shown to be consistent with a kappa value of two [123]. 
The Qn resonances also have subsidiary resonances [124, 
125]. For a particular n value, the Qn resonance and its 
subsidiary have been attributed to the Doppler shifts of 
two waves, one propagating as a forward mode and the 
other as a backward mode, near but on both sides of the 
un-shifted Qn frequency [121]. The first member of each 
sequence is shown in Figure 4 [67]. This also shows one of 
the D 1 subsidiary resonances (D 1 +); the hybrid relationship 
proposed to explain it [ 1 1 9] ; and a resonance designated as 
“F” because it is often observed to float from the zero time- 
delay (or zero apparent range) baseline, and the approximate 
expression for its frequency [126]. 
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Figure 4. An ISIS-2 digital ionogram illustrating the Dl, D1+, Q2, and F resonances (in addition to the 
resonances that can be explained by the oblique-echo mechanism characterized by Figure 3) [67], 


The sequence nature of the Dn and Qn resonances 
is not apparent from a single ionogram, such as the one 
reproduced in F igure 4. However, as the plasma parameter, 
fpej fee ’ increases, higher members of the sequences 
appear. A presentation displaying the observed resonant 
frequencies from many ionograms, normalized by the 
observed f ce versus the observed f pe jf ce illustrates 
the dominant control of the plasma parameter f pe j f ce in 
determining the frequency of these resonances [108, 118, 
120, 125]. These resonances have proven to be of great 
value in the spectral classification of sounder-stimulated 
plasma resonances for the accurate determination of N e 
, particularly in the magnetosphere, where clear X and Z 
traces are often not present [68]. While the Dn-sequence 
portion of this spectral classification is well explained 
by the cylindrical-electromagnetic plasma oscillations 
mechanism, some major issues remain. For example, these 
issues include deriving the fundamental D 1 frequency from 
first principles (presently it is based on an empirical fit to 
the observations), and incorporating propagation aspects 
likely needed to explain the observed long time durations 
of the Dn resonances and the frequency splitting of the D 1 
resonance observed at mid latitudes [108, 127]. The issues 
also include explaining the abrupt disappearance of the 
3 f ce resonance as the D2 resonance becomes prominent 
[125,128] (see also Section 6), and investigating the role of 
sounder-stimulated ion motions leading to ion signatures on 
Dn resonances [128] (such ion signatures will be discussed in 
Sections 5 and 6). These issues, plus many other features of 
the Dn and Qn resonances, such as their excitation by natural 
process and (in the case of Dn) by rocket-borne electron 
guns, their use in determining N e , and their relationship 
(in the case of Dn) to sounder-accelerated electrons, are 
discussed in Sections 1 and 3 of [111], Sections 3 of [68] 
and [113], and the “Discussion and Summary” of [129]. 
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The revolutionary concept of interpreting the plasma 
resonances as oblique echoes, introduced by McAfee [100] 
after building on the workofCalvert [ 1 30] in interpreting the 
Z' echo (see, e.g., Figure la), was confirmed by innovative 
rocket experiments, as described in the next section. 

4. Instrumental Innovation 

Having had the privilege of worldng with Willi 
Hough on the ionosphere program, using the massive C3 
ionosonde, during the first winter-over operation in 1957 
at the Amundsen-Scott 1GY South-Pole Station, under the 
leadership of Paul Siple [131], the author fully appreciates 
the degree of instrumental innovation involved in placing 
a swept-frequency ionospheric sounder on a satellite in 
1962 [4, 55, 132, 133]. At the time, many thought that the 
Canadian-built Alouette-1 would at best have the then- 
typical satellite lifetime of a few months. The pessimism 
that the sounder would work was so strong that NASA had 
made no plans to use data from it. However, thanks to a 
Robust Design approach, the sounder provided a wealth 
of data for 10 years, even surviving higher-than-expected 
radiation from the unanticipated high-altitude hydrogen- 
bomb test over Johnston Island in the Pacific Ocean, seven 
weeks before the Alouette-1 launch [134, 135]. 

Several rocket payloads in the early 1970s were 
designed specifically to investigate the plasma resonances 
detected by ionospheric topside sounders. Rather than 
retrieving only the received pulse envelope, as in the case 
of the Alouette and ISIS topside sounders, these rocket- 
borne sounders telemetered the sounder-receiver waveform 
to the ground. Fourier analysis of this waveform resulted 
in unprecedented frequency accuracy. In one experiment, 
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Figure 5a. The Fourier spectrum of the 3 f ce resonance waveform 
of the signal observed at 3.096 MHz when the rocket-borne sounder 
was at 480 hn altitude, where f pe = 4.20 MHz [7], 

the 3 f ce resonance was shown to be the result of a beating 
of two waves separated in frequency by about 1 kHz, as 
illustrated in Figure 5a [7]. These remarkable experimental- 
rocket results suggested that the 3 f ce resonance was the 
result of oblique echoes of electrostatic waves, i.e., the 
3 f ce resonance was the result of a mechanism similar to 
the one proposed earlier for the f and f uh resonances 
[100, 101, 136]. In approximately the same timeframeas this 
rocket experiment, ray-tracing calculations demonstrated 
that an Alouette-2 3 f ce resonance could be explained 
by the oblique-echo mechanism [102]. In another rocket 
experiment, the frequency variation of the received f uh 
signal with delay time was measured, and found to agree with 
theoretical predictions based on oblique echoes of slowly 
propagating sounder-generated waves (see Figure 5b). 
Previous satellite ionospheric topside sounders could not 
observe this frequency variation with delay within the 
receiver bandwidth, because of onboard signal detection 
[11]- 

The highly successful bistatic rocket sounders 
OEDIPUS A and C obtained magnetic-field-aligned 
separations of the order of 1 km in the auroral topside 
ionosphere with the sounder transmitter on one payload and 
the sounder receiver on the other. OEDIPUS-C provided 
the first in-situ demonstration of Faraday rotation in space 
[30]. It also provided one of the first controlled quantitative 
confirmations of incoherent radiation theory in space, by 
detecting slow Z-mode radiation from sounder-accelerated 
electrons [137]. The OEDIPUS-C measurements of field- 
aligned irregularities, responsible for wave ducting of 
sounder generated X- and Z-mode waves, showed they 
consisted of N e depletions of up to 21% with dimensions 
across B of a few kilometers [138]. 


Figure 5b. The echo delay as a function of frequency 
for the f u h resonance when f u ^ was slightly less than 
2 fee ■ different symbols indicate different combi- 
nations of transmitter power and receiver gain during 
different modes of operation triggered by the large 
amplitude of the f u ^ resonance, and the solid lines 
are theoretical predictions based on different electron 
temperatures [11]. 

The above rocket experiments introduced innovations 
in sounder design that made use of high data rates on short- 
duration missions dedicated to radio sounding. The design 
considerations for the sounder known as the Radio Plasma 
Imager (RPI) for the magnetospheric IMAGE satellite were 
quite different. IMAGE was a remote-sensing observatory 
with many instruments, and was far from being a dedicated 
sounding mission. The Radio Plasma Imager had to 
overcome programmatic challenges as well as instrumental 
challenges [139]. The instrumental challenges were many. 
Sounding was to be performed over great distances involving 
reflections from low-density plasmas, in spite of power, 
weight, and data-rate restrictions [66, 140]. The solution 
was to use 500 m tip-to-tip dipole antennas, transmit coded 
pulses that enabled signal processing to increase the signal- 
to-noise ratio so as to obtain long-range echoes with only 
10W radiated power, and to have a variety of programs 
and schedules to optimize the science return within the data 
restrictions [33]. The scientific results have been profound, 
as illustrated in a number of Radio Plasma Imager review 
papers [34, 141-143], reviews based on IMAGE/RPI and 
Cluster/WHISPER [144-147], and in other sections of 
this review. 

Similarly, the design considerations forthe WHISPER 
relaxation sounders on the magnetospheric Cluster satellites 
were also quite different. Here, the challenge was to deliver 
simultaneous N e measurements on board the four spacecraft 
of the constellation in order to calculate instantaneous 
N e spatial gradient vectors. The variability of plasma 
conditions nearboundaries demanded a duration of each N e 
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Figure 6. Six consecutive ISIS-1 
ionograms when the sounder 
was operating in the G mode 
(alternating between fixed-fre- 
quency and normal fixed/swept- 
frequency ionograms), illustrat- 
ing the antenna spin modulation 
of the Dl* resonance, the 2 fjj 
resonance, and its ion spur. The 
optimal spin phase for reso- 
nance excitation, the calculated 
value for the Dl resonance, and 
the calculated value for the Dl* 
resonance are designated at the 
top of each ionogram by * Dl, 
and +, respectively. 


determination significantly shorter than had been obtained 
on board previous magnetospheric sounders, like GEOS or 
ISEE. This requirement was met by implementing a fast 
Fourier transform, instead of a swept-frequency analyzer, 
in the receiver and analysis chain. The time resolution thus 
obtained ( — 1.5 s) allowed quasi-synchronization of the 
four independent measurements, via on-board UT clocks 
[35, 122, 148, 149], 

5. Sounder-Stimulated Plasma 
Resonances: Importance of the 
Antenna’s Orientation 

There are several factors that can determine the 
occurrence of, and the characteristics of, sounder-stimulated 
plasma resonances, as have been discussed in a number of 
reviewpapers[61, 102, 103, 150, 151]. Lockwood [152] was 
the first to show the importance of the antenna’s orientation 
on the sounder excitation of plasma resonances. He showed 
that nf ce resonances with high n were only observed by 
the Alouette- 1 topside sounder when the relevant sounder 
antenna was nearly parallel to B. The relevant antenna 
was the one being used for transmission and reception. 
(In order to cover the wide frequency range, Alouette- 1 
employed crossed orthogonal dipoles of different lengths 
in the spin plane, with the longer dipole being used for 
transmission and reception at the lower frequencies, and 
the shorter dipole used for transmission and reception at 
the higher frequencies. The Alouette-2, ISIS-1, and ISIS-2 
sounders operated in a similar manner [4].) While the 
effective radiated power does not appear to be a significant 
factor in determining the duration of the nf ce resonances, 
as indicated by a comparison of selected Alouette- 1 and 
Alouette-2 observations (see Section 5 of [106]), plasma 
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conditions do play an important role. The following three 
studies indicated that for large f pe j f ce , plasma conditions 
may be the dominant factor: 

1 . An investigation of the time durations of Alouette-2 
nf ce resonances as a function of altitude near the dip 
equator found the greatest durations for resonances 
with n > 5 in the relatively dense plasma near perigee 
(see Figures 1 and 8b of [106]), where large values of 
fpe /fee are expected. 

2. An inspection of more than 200 Alouette-2 ionograms, 
selected so as to provide a wide range of f pe J f ce values 
within a narrow range of f ce values, found few high- 
order nf ce resonances excited when f pe /f ce < 4 . but 
found the number of high-order harmonics to increase 
linearly with increasing f pe j f ce (see Figures 2, 21, 
and 22 of [125]). 

3 . Adisplay of ionograms from two ISIS-1 passes, selected 
to obtain a range of f pe j f ce values from approximately 
one to 10, again while maintaining a narrow range 
of f ce values, revealed a striking decrease in the 
duration and frequency width of the nf ce resonances 
when nf ce > f uh / f ce (most apparent for n > 2 ). The 
greatest nf ce resonant duration corresponded to the 
resonance with frequency closest to, but less than, / 
[153]. Thus, these studies suggested that a particular 
antenna orientation is not required for the excitation 
of high-order nf ce resonances when f pe jf ce is large 
( > 4 ). These results were consistent with the antenna- 
orientation restrictions for the detection of high-order 
nf ce resonances found by Lockwood [152], since the 
Alouette- 1 data set he used only sampled f pe jf ce 
values less than about four [125]. 
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The sounder’s antenna orientation plays an important 
role in determining the characteristics of the lower-order 
nf ce resonances, but that role appears to be more complex 
than in the case of the high-order nf ce resonances when 
fpef fee * s sma ll (< 4 ). An investigation of Explorer 20 
fixed-frequency observations indicated that the 3 f ce 
resonance revealed a peak in duration with a pronounced 
fringe pattern when the antenna was perpendicular to B, 
but even stronger and longer-duration signals (without 
fringes) when it was parallel to B (see Figure lb of [6]). 
The striking fringe pattern observed when the antenna 
was perpendicular to B is suggestive of two waves beating 
together, as observed for the 3 f ce resonance in a rocket 
experiment [7] as the antenna orientation relative to B varied 
approximately ±20° around the perpendicular direction. 
The rocket observations were explained in terms of oblique 
echoes using analytic ray tracing [154]. This analytic ray- 
tracing approach was applied to Alouette-1 low-order nf ce 
resonance observations, where it was found that near the 
magnetic equator, the maximum resonant duration would 
be expected when the antenna was nearly perpendicular 
to B (see Equation 31 of [104]). Low-latitude Alouette-1 
observations indicated that strong plasma resonances at / 

, f uh , 2 f ce , and 3 f ce could be observed using the long 
antenna when it was far from parallel to B, since high-order 
nf ce resonances were observed a few seconds later on the 
short antenna [1 55] (an assumption supportedby the plasma 
conditions, which corresponded to f pe j f ce < 4 ). ISIS-1 
fixed-frequencyobservationsofthe f resonance revealed 
strong maxima when the antenna was nearly perpendicular 
to B that were separated by a deep minimum when it was 
perpendic ular to B. These observations have been explained 
by including the finite size of the antenna in the oblique- 
echo interpretation [ 1 56] . In general, it appears that plasma 
resonances at / , f uh , 2 f ce , and 3 f ce can be stimulated 
over a wide range of antenna orientations. However, the 
signal strengths, time durations, and fringe patterns change 
in spectacular ways when the angle between the antenna 
and B is near either 0° or 90°. 

Figure 6 presents six consecutive ISIS-1 ionograms, 
illustrating the enhancement of the 2 f ce resonance (the 
notation 2 f H is used in the figure) when the sounder antenna 
was in an optimal orientation. During this sequence, the 
ISIS- 1 topside sounder alternated between fixed-frequency 
ionograms (Figures 6a, 6c, and 6e) and normal fixed/ 
swept-frequency ionograms (expanded portions shown 
in Figures 6b, 6d, and 6f). From the later ionograms, it 
was clear that f pe j f ce < 3 , and that the 2 f ce resonant 
frequency was significantly greater than 1.0 MHz, the 
frequency corresponding to the fixed-frequency operations. 
Note that even though the 2 f ce resonant frequency was 
greater than 1.0 MHz, the 2 f ce resonance was observed 
on the fixed-frequency ionograms (Figures 6a, 6c, and 
6e) for a short duration with an intense (red) signal every 
time the sounder antenna was in an optimum orientation 
(these orientations are indicated by the * marks at the top 
of the ionograms). These short-duration intense signals 
resulted from a broadening of the bandwidth of the 2 f ce 
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resonance during these favorable antenna orientations. 
(These signals shifted in phase from one fixed-frequency 
ionogram to the next fixed- frequency ionogram because the 
time interval between these ionograms was near but not an 
exact multiple of the antenna-spin period, the former being 
57.9 s and the latter being about 20 s.) This broadening (to 
include 1 .0 MHz) was seen on only one of the normal fixed/ 
swept-frequency ionograms, namely, Figure 6d, the only 
ionogram where the sounder antenna was in an optimal 
orientation as the swept frequency was near 1 .0 MHz (as 
indicated by the * in Figure 6d). In addition to the larger 
bandwidth, the intense (red) portion of the 2 f ce resonance 
was observed to have a longer duration on this ionogram 
than on Figures 6b and 6d. These characteristics suggested 
that the optimal antenna orientations, flagged by the * 
marks, corresponded to the antenna being nearly parallel to 
B, based on the Explorer-20 fixed-frequency observations 
where the strongest nf ce resonances with no fringe patterns 
were observed under these conditions [6, 150]. 

There were two other striking features to notice in the 
expanded portion of the fixed/swept-frequency ionogram 
shown in Figure 6d, where this optimal antenna orientation 
occurred when the swept frequency was near 1 .0 MHz. F irst, 
the 2 f ce resonance had an ion spur near 500 km apparent 
range that extended far enough on the low-frequency side 
of the resonance so as to reach the 1.0 MHz frequency 
marker. Second, there was a long-duration D1+ resonance 
near 1.0 MHz extending to nearly 500 km apparent range. 
Both of these features were only apparent on the fixed- 
frequency ionograms (Figures 6a, 6c, and 6e) under the * 
marks, indicating optimal antenna orientation, and they were 
less pronounced during the swept-frequency operation on 
Figures 6b and 6f, where the optimal antenna orientation 
(* marks) did not occur near the excitation of the 2 f ce 
resonance or the calculatedpositions ofthe D 1+ resonances 
(indicated by the arrows below the + signs near 1 .0 MHz). 
No theoretical explanation has been proposed to explain 
this dependence of the D1+ resonance stimulation on the 
antenna’s orientation. However, this dependence was likely 
the reason the subsidiary Dn resonances were not always 
present with their corresponding main resonance on the 
same ionogram. The example in Figure 4, where the D 1 and 
D1+ resonances were clearly seen on the same ionogram, 
is not typical. 

The ion spur illustrated in Figure 6 (also known as a 
proton spur, because it appears after a time delay equal to 
the proton gyro period) clearly demonstrated that positive 
ion motions cannot be ignored if such sounder-stimulated 
phenomena are to be explained. These spurs were also 
apparent in the ionogram ofF igure 4 near 500 and 1 ,000 km 
apparent range on the f pe resonance (labeled “N”), and on 
the 2 f ce resonance (labeled “2”), where the spur extended 
nearly to the 1 . 0 MHz frequency marker at an apparent range 
slightly below the 500 km apparent-range label. Evidence for 
sounder-stimulated ion motions was first identified as a series 
of spurs on the sounder-generated plasma resonances on 
Alouette- 1 topside ionograms [157]. These spurs have been 
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the subject of additional investigations using later topside 
sounders with better frequency resolution that identified 
specific resonances containing ion spurs [128, 158, 159]. 
Their explanation may be similar to the interpretation of the 
proton-cyclotron echoes observed on topside ionograms. 
These interpretations involve either proton bunching 
due to the RF -pulse-induced negative sounder-antenna 
potential [160, 161], or selective proton acceleration by 
the RF sounder pulse depending on the RF phase [162, 
163]. Recent IMAGE/RPI observations at various altitudes 
within the plasmasphere have revealed intense W-mode 
proton-cyclotron echoes that were attributed to proton 
excitation as a transient event near the beginning of each 
RF sounder pulse. They were believed to be the most 
efficient at frequencies near and below the proton plasma 
frequency [57]. These Radio Plasma Imager observations 
also revealed two other forms of proton-cyclotron echoes 
that were similar to echoes observed by the Alouette/ISIS 
topside sounders. The first of these involved striking echoes 



Figure 7. Portions of three consecutive ISIS-1 ionograms 
illustrating the changes in sounder-stimulated plasma 
phenomena as the ambient plasma conditions changed from 
f pel fee < 5 in ( a ) t0 f pel fee > 8 in ( c ')■ Spectacular 
enhancements occurred when f pe / f ce ~ 6 near 13 s after 
the frame sync, corresponding to a swept-frequency value 
of about 4.0 MHz in (b). (Recorded at QUI on 1971_006 
with ISIS-1 at -11.9 GMLONG and 1914 GMLTM; 0037:36, 
0038:06, and 0038:05 UT, 13.4, 11.5, 9.6 GMLAT, 704.5, 
688.1, 673.3 km, in (a), (b), and (c), respectively: scaled 
f pej fee values: 3.03/0.676. 3.97/0.656 . and 5.45/0.642 
in (a), (b), and (c), respectively.) 


at frequencies just above f ce with a delay-time dependence 
on sounder frequency, suggesting warm-plasmapropagation 
following a sounder-perturbed proton distribution possibly 
similar to that proposed for the ISIS observations [162]. 
The second involved proton-cyclotron echoes associated 
with the cold-plasma Z mode [57]. 

6. Sounder-Stimulated Plasma 
Resonances: Importance of 
Special Plasma Conditions 

Sounder-stimulated plasma phenomena have been 
observed to be particularly sensitive to the value of 
the plasma parameter f pe /f ce ■ There was a stepping 
through higher members of the Dn resonances as f pe j f ce 
increases [108, 118, 125]; an instability signature on the 
nf ce resonances when nf ce < /„/, [153]; an enhancement 
of the proton spurs when nf ce ~ f pe [128]; a strong 
diffuse resonance between / and f uh (known as DNT) 
that changed from a floating to a non-floating resonance 
as the plasma conditions changed from / pe // ce < 1 to 
fpe / fee > 1 [125]; an abrupt disappearance of the 3 f ce 
resonance when f pe j f ce a 4 (see Figure 9 of [128] and 
Figure 14 of [125]); and enhanced scattering of sounder- 
generated Z-mode waves when f pe l f ce « n and n > 4 , 
suggesting either sounder-enhanced or sounder-generated 
field-aligned irregularity [164]. Figure 7b illustrates this 
last phenomenon. Here, enhanced proton spurs and Z-mode 
scattered signals were observed when f pe / f ce « 6 . During 
the three-ionogram sequence of Figure 7, f ce decreased 
slightly from F igure 7a to 7c, and thus the frequencies of the 
nf ce resonances varied only slightly from one ionogram to 
the next. Some of the more -prominent nf ce resonances in 
Figure 7c are labeled at the top of the ionogram. However, 
the / and f uh resonances (coalescing into the long- 
duration signal between the Z and X reflection traces) 
increased considerably in frequency from one ionogram 
to the next. Only in Figure 7b, where the indicated signal 
enhancements were observed, was the ratio f pe /f ce near 
an integer value. 

It has been shown that the enhanced signals attributed 
to wave scattering in F igure 7b were not due to more-efficient 
scattering when f pe j f ce « n [165]. This result supports 
the argument that the wave scattering is from field-aligned 
irregularity stimulated, or enhanced, on a short time scale 
by the sounder [164]. These observations stimulated two 
theoretical studies: 

1 . An investigation of large -amplitude electron oscillations 
in cylindrical plasma structures aligned along B, i.e., with 
structures that could represent field-aligned irregularity, 
where enhanced resonance conditions were found when 
fpe / f ce ~ n and where it was suggested that future 
work include dissipation and driver mechanisms [ 1 66]. 
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Figure 8. An IMAGE/RPI plasmagram (color-coded signal amplitude, with red most intense, displayed as virtual 
range from 0 to 7.5R E on a linear scale as a function of sounding frequency from 400 to 1000 kHz on a logarithmic 
scale) with multiple Z- and X-mode epsilons, each containing traces due to multiple-hop hemisphere-to-hemisphere 
echoes. Echoes within the local hemisphere, i.e., the hemisphere containing the IMAGE satellite, are labeled “S” for 
short echo path, and echoes from the conjugate hemisphere are labeled “L "for long echo path. The outside labels on 
the bottom scale designate the Z- and X-mode cutoffs. The inner three labels on the bottom scale, and the three arrows 
at the top, identify the model value for 2 f ce and the calculated f and f u ^ values consistent with the observed 
wave cutoffs. (1850:28 UT on 31 August 2002; 11.1 MLAT, 0228 MET, l.5R E , L =1.57 .) 


2. A four-wave coupling process, where a sounder- 
generated Z-mode wave acts as a pump wave at 
frequency f 0 , which decays into two electron Bernstein 
waves and a purely growing meter-scale field-aligned 
irregularity when this Z-mode wave propagates to the 
O-mode reflection layer, where f a ~ f pe ~ nf ce [167]. 
In this latter study, the f « nf ce condition would thus 
be achieved at some distance from the satellite, rather 
than at the satellite’s location as observed. 

A new resonance, stimulated by the Radio Plasma 
Imager on the IMAGE satellite, has been observed at 
a frequency about 15% above f ce when the spacecraft 
is within the plasmasphere but at altitudes above about 
7,000 km [57] . It was suggested that the lack of observation 
of this resonance by the Alouette-2 and ISIS-1 topside 
sounders may have been due to their relatively low apogees 
of 3,000 and 3,500 km, respectively [57]. Recall that in 
Section 3 , it was pointed out that the there is both theoretical 
and observational evidence that the f ce resonance itself is 
also expected to be observed in a stronger form in a rare 
rather than a dense plasma. 

The plasma conditions encountered by the URAP 
relaxation sounder on Ulysses in J up iter ’s Io plasma torus led 
to a spectrum of sounder-stimulated plasma resonances quite 
different from those encountered in the Earth’s ionosphere 
or magnetosphere. The resonances were first interpreted in 
terms of the Dn resonances and the / resonance, with no 
evidence for resonances at nf ce or f uh , to deduce both |fi| 
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and N e along the spacecraft(s/c)trajectory [27, 168]. Later, 
they were interpreted as Doppler-shifted nf ce resonances, 
due to the large (> 100 km/s) spacecraft velocity relative 
to the co-rotating Jovian magnetospheric plasma, and were 
used to deduce both N e and T e along the spacecraft’s 
trajectory [28]. These different interpretations generated 
some controversy [169-171]. 

7. Nearly Lossless Propagation 
within N e Wave Ducts 

The first satellite evidence of efficient propagation 
within N e wave ducts was in the classic paper by Muldrew 
[172], published one year after the Alouette-1 launch. 
There have been many investigations of this non-vertical 
propagation [86] that include the concept of ducting within 
field-aligned irregularity in the whistler (W), Z, O, or X 
modes [83, 87, 173-176]. These ducted signals can take 
the shape of spectacular signals appearing as floating 
epsilons on ionospheric topside ionograms [177, 178] 
and on magnetospheric plasmagrams [83, 176, 179], as 
illustrated in Figure 8. 

The intensities of the epsilon signatures, which are 
caused by multiple hemisphere-to-hemisphere echoes, 
suggest nearly lossless propagation within wave ducts. Such 
propagation has been well described for the W mode, as 
well as for the other modes (see, e. g. , [ 1 73 , 1 74] ) . Conj ugate 
O-mode epsilons are not apparent in this figure, but when 
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they are observed, they provide an independent confirmation 
of field-aligned N e profiles deduced from Radio Plasma 
Imager plasmagrams. This confirmation has indicated that 
what appears as an O-mode magnetospheric conj ugate echo 
is the result of O-Z-O-wave-mode coupling [179], rather 
than direct O-mode ducting as proposed to explain the 
Alouette-1 conjugate echoes observed in the ionosphere 
[172]. The subject of O-Z mode coupling is also of great 
interest in ground-based wave-injection experiments (see, 
e.g., [ 180 , 181 ]). 

It is to be emphasized that the combination of echoes 
from the conjugate hemisphere with echoes from the local 
hemisphere enabled nearly instantaneous magnetic-field- 
aligned N e profiles to be produced from one hemisphere to 
the other (with only the low-density region near the equator, 
where N e is less than the value at the IMAGE satellite, 
being filled in by interpolation). As IMAGE travels along its 
orbit, consecutive hemisphere-to-hemisphere field-aligned 
profiles can be obtained for different radial distances, 

i.e., different L values. Such consecutive Radio Plasma 
Imager field-aligned N e profiles have been combined 
to produce orbit-plane plasmaspheric N e contours for a 
number of IMAGE orbits before, during, and after the large 
magnetic storm of March 31, 2001, in order to describe 
the plasmaspheric mass loss and refilling processes [182]. 

8. Summary 

Space -borne radio sounding is a highly reliable 
technique, compatible with other instruments, for accurately 
measuring the local |fi| , and both local and remote N e . 
Sounder observations have provided a wealth of information 
on the terrestrial topside ionosphere and magnetosphere, 
and on basic processes of wave generation and propagation. 
In spite of these great achievements - which have also 
been obtained in extraterrestrial plasmas - additional 
improvements to the International Reference Ionosphere 
for the topside ionosphere are required to meet the needs of 
present-day technological systems. Fundamental questions, 
such as the following, remain to be answered concerning 
many sounder-stimulated plasma phenomena: 

1 . What processes, in addition to the oblique-echo model, 
are needed to explain satellite observations of the 2 f ce 
resonances? 

2. What is the cause of the abrupt disappearance of the 
3 f ce resonance when f pe / fee ® 4 ? 

3. Can the oblique-echo model explain the long durations 
observed for the magnetospheric nf ce resonances? 

4. What is the origin of the plasmaspheric resonance 
observed about 1 5% higher in frequency than f ce , and 
the proton-cyclotron echoes observed both in a similar 
frequency region and near f 7 ? 

5. Will mechanisms applied to ionospheric and 
magnetospheric proton echoes be able to explain the 
proton spurs observed on sounder-stimulated plasma 
resonances, including their dependence on f pe jf ce 
and antenna orientation? 
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6. Can the fundamental D1 frequency be derived from 
first principles, can propagation aspects be incorporated 
into the cylindrical-electromagnetic plasma oscillations 
mechanism so as to explain the observed durations of 
the Dn resonances, and can the antenna spin-modulation 
dependence of the subsidiary Dn resonances be 
explained? 

7. What causes the strong DNT resonance observed 
between f pe and f uh when f pe jf ce is near but less 
than one? 

8. How does a space-borne sounder enhance, or 
generate, field-aligned irregularity when f pe J f ce « n 
and n > 4 ? 
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